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DETECTION  OF  ELEMENTARY  PROCESSES  WITHIN  POROUS  ELECTRODES 


S.  S?pak,  T.  Kdtdn 

Naval  Ocean  Systems  Center,  San  Diego,  CA  92152 
Lockheed  Missiles  & Space  Co.,  Inc. 

Lockheed  Palo  Alto  Research  Lab,  Palo  Alto,  CA  94304 

ABSTRACT 

Problems  in  electrode  modeling  are  reduced  if  parti- 
cipating processes  are  recognized  and  experimentally 
defined.  Difficulties  in  treating  complex  interplay 
between  elementary  processes  within  confines  of  a changing 
porous  structure  are  minimized  through  application  of 
simple  concepts,  e.g.,  penetration  depth  and  a structural 
unit  characteristic  of  a porous  electrode.  A systematic 
approach  to  the  determiriation  of  relevant  processes  and  to 
specification  of  electrode  structure  for  a prescribed  per- 
formance is  examined.  The  analysis  is  applied  to  the 
Ag/AgCt  system  in  1 N KCt  operating  at  23°C. 

In  work  concerned  with  electrodes  for  electrochemical  systems, 
particularly  in  the  field  of  energy  conversion  and  storage,  two  deci- 
sive design  factors  have  emerged,  reaction  intensity  profile  and 
thickness  of  the  porous  electrode.  This  generalization  applies 
equally  to  fuel  cells,  to  various  hybrid  systems  as  well  as  to  conven- 
tional batteries.  This  seems  to  imply  that  a rational  electrode  design 
should  proceed  via  an  electrode  modeling  approach. 

The  mathematics  of  electrode  modeling  has  advanced  to  a high 
degree  of  sophistication.  Today's  models,  chiefly  due  to  efforts  of 
Bennion  and  Newman  and  their  co-workers (1 ) , are  applicable  to  complex 
situations.  These  models  are  capable  of  handling  such  factors  as 
changes  in  surface  morphology,  reaction  kinetics  and  the  transport 
properties.  These  factors  however,  must  be  determined  experimentally 
so  that  correct  computer  inputs  can  be  provided  for  the  various  sys- 
tems under  consideration.  The  present  status  may  be  summarized  as 
follows:  the  mathematics  of  electrode  modeling  are  well  developed  but 

the  realistic,  experimental,  inputs  required  for  meaningful  modeling 
are  lagging  behind. 

It  is  the  purpose  of  this  coimiuni cat ion  to  examine  methods  by 
which  elementary  processes  operating  within  the  porous  structure  can  be 
detected  and  identified.  The  emphasis  is  on  the  use  of  a concept  of 
penetration  depth  and  the  examination  of  surface  morphology.  The 
analysis  is  applied  to  the  Ag/AgCt  system  in  1 N KCt  operating  at  room 
temperature. 
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ELECTRODE  ANALOGS  AND  SIMULATION  OF  OPERATION 

Studies  of  the  dynamics  of  porous  electrode  operation  suggest 
that  structural  analogs  can  be  constructed  to  provide  useful  tools  in 
guiding  the  development  of  practical  battery  electrodes. (2 ,3)  To  il- 
lustrate, consider  Fig.  1 where  a set  of  elementary  pr">cesses  usually 
associated  with  an  electrochemical  reaction.  Fig.  lb,  are  taking  place 
within  the  confinements  of  a pore.  Fig.  la.  If  we  accept  the  premise 
that  the  characteristics  of  a porous  electrode  can  be  represented  by  a 
set  of  events  on  planar  electrodes  which  operate  subject  to  local  con- 
straints on  the  reaction  path,  then  the  search  for  correct  inputs  for 
the  electrode  modeling  may  be  accelerated. 

Consider  Fig.  2 to  relate  the  experimentally  simulated  reaction 
intensity  profile  to  that  associated  with  the  operation  of  a real 
porous  electrode.  In  particular.  Fig.  2a  represents  a real  electrode 
structure.  For  the  purpose  of  construction  of  an  analog,  we  select  a 
volume  element  sufficiently  la'~ge  compared  to  the  individual  grains, 
but  small  enough  compared  to  the  penetration  depth.  The  electric  cir- 
cuit analog  is  shown  in  Fig.  2b  while  the  experimental  arrangements, 
or  the  structural  analogs,  are  shown  in  Figs.  2c  and  2d.  Specifically, 
Fig.  2c  shows  a small  section  of  an  electrode  structure  generated  by  a 
loose  packing  of  metallic  spheres  while  Fig.  2d  illustrates  the  prin- 
ciple and  construction  of  a segmented  electrode. 

The  concept  of  an  electric  circuit  analog  is  introduced  here  to 
provide  a basis  for  an  elementary  treatment  of  porous  electrodes. 

Thus,  it  is  assumed  that  in  real  porous  electrodes,  the  position  de- 
pendent current  i{x)  and  potential  u(x)  in  the  electrolyte  phase  are 
related  by  Eqs . (1)  and  (2) 
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= 0 


3i  u 
8x  Z 


= 0 


(1) 

(2) 


The  solution  of  Eqs.  (1)  and  (2)  exhibits  an  exponential  behavior  with 
a curvature  which  is  dependent  on  the  numerical  value  of  = R/Z. 

For  a sufficiently  thick  electrode,  the  solution  is 

y(x)  ^ A exp(-Kx)  (3) 

where  y(x)  stands  for  either  i(x)  or  u(x)  and  A is  determined  from  the 
applicable  boundary  condition.  However,  because  j = -3i/3x  and,  for  a 
conductive  matrix,  u(x)  = n(x)  const.,  ft  follows  that  y(x)  yields 
the  reaction  intensity  profile,  j(x),  as  well  as  the  distribution  of 
the  driving  force,  the  overpotential , n(x). 

The  character! Stic  number,  ic,  governing  the  distribution  functions 
within  the  porous  structure  has  been  converted  to  a useful  engineering 


number.  It  has  been  used  in  tne  design  chemical  reactors  since 
1939  and  was  defined  simply  as  the  ratio  ot  transport  resistance  to 
the  reaction  impedance . (4)  It  was  in  tliis  form  that  Bro  and  Kang(5) 
applied  it  to  describe  the  behavior  of  sintered  cadmium  electrode. 

In  electrochemical  engineering  its  significance  was  clearly  expressed 
by  Wagner(6)  and  it  gained  acceptance  in  the  1960's  in  connection  with 
the  design  of  fuel  cell  electrodes. 

A selected  set  of  the  characteristic  number,  together  with 
respective  author's  formulation  and  nomenclature,  are  presented  in 
Table  I. 


TABLE  I. 

Selected  Examples  of  Formulation  of  the  Characteristic  Number 


Author  Formulation  Nomenclature 


Bro  and  Kang  (5) 

ionic  impedance  _ „? 
faradaic  impedance 

electrochemical 
Thiele  modulus 

Wagner  (6) 

^ - oRT 

^C  a zFJ 
c 

polarization  param- 
eters 

Nanis  (7) 

M9  _ 1 9n  1 

^ L 3j  • LTW 

throwing  power 

Grens  and  Tobias  (8) 

^ ^ io 

nFDkCk° 

transfer  current 
distribution  param- 
eter 

Gidespow  and 
Baker  (9) 

1 

n io  nF  aL^ 

^ RT  K 

dimensionless 
exchange  current 

Alkire  and 
Mi  rare fi  (10) 

2nFfc2  i^ 

' V "o" 

dimensionless 
reaction  velocity 
parameter 

Winsel  (11) 

^0 

effective  depth  of 
penetration 

Physical  meaning  of  the  characteristic  number.  The  characteris- 
tic number  can  be  formulated  either  as  a dimensionless  number,  in  a 
form  indicating  the  curvature  of  the  distribution  function,  or  as  a 
measure  of  penetration  of  the  reaction  zone  into  the  electrode 
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structure.  It  also  can  be  viewed  as  a product  of  quantities  associated 
with  the  electrode  reaction,  inclusive  of  transport  and  properties  of 
electrode  structure.  This  is  apparent  in  Nanis' (7)  formulation. 


To  illustrate  the  kinetic  nature  of  the  penetration  depth  and  to 
relate  it  to  the  electrode  structure,  consider  Fig.  2d  and  apply 
Eq.  (1)  to  an  element  dx,  in  a form 


ai (x)  = -o  adu/dx  (4) 

where  a is  the  spacing  between  the  electrode  segment  and  the  glass 
cover,  and  o is  the  electrolyte  conductivity.  Eq.  (4)  together  with 
an  expression  for  charge  conservation 

d(ai)/dx  + j = 0 (5) 


yields,  upon  substitution  of  Eq.  (5)  into  Eq.  (4),  a differential 
equation 


dx 


- j = 0 


(6) 


Assuming  for  simplicity  a constant  a and  a,  and  linearizing  the  Tafel- 
like  j(n)  relation,  we  obtain  the  familiar  equation 


d^n 

dx^ 


RT 

oaF 


n = 0 


(7) 


and  therefore  also  an  expression  for  the  characteristic  length  which 
can  be  written  in  the  form  of  a product  of  the  electrode  reaction 
parameters  and  the  geometry  of  segmented  electrode,  namely 


a 


(8) 


It  is  easy  to  see,  cf . , Fig.  2d,  that  a = V/S,  where  V is  the 
electrolyte  volume  and  S the  interphase  area.  For  an  arbitrary  elec- 
trode structure,  a is  replaced  by  ()>  = <V>/<S>,  cf.  , Fig.  2c.  The 
structural  parameter,  ip,  will  remain  constant  as  long  as  the  structural 
unit  is  maintained. (12) 


The  effect  of  a change  in  the  reaction  path  on  the  electrode  be- 
havior can  now  be  examined  in  terms  of  the  characteristic  length.  To 
account  for  any  additional  step,  the  so-called  approach  step(13),  we 
may  write  an  expression 


/I 


/FtRT 


(9) 


The  less  familiar  symbol,  i , is  the  phenomenological  coefficient  which 

relates  the  reaction  velocity,  v,  to  the  driving  force,  A.  To  account 
for  the  set  of  contributing  processes,  we  proceed  as  follows:  the 

*\j 

electrochemical  affinity.  A,  of  the  overall  process  is  taken  to  be 
equal  to  the  sum  of  affinities  of  the  individual  elementary  processes, 
i .e. , 


A = E (k  = 1.  2,  3,....) 


(10) 


By  Eq.  (10),  only  a set  of  consecutive  reactions  is  allowed.  For 
a stationary  state:  v^  " ^2  ~ ~ ^k  “ because  of  the  rela- 

tionship V|^  = t|^A|^,  we  have  immediately 

£-1  = 2 (”) 
k 


Substitution  of  Eq.  (11)  into  Eq.  (9),  yields 


j = -I 


axp(-  - 


1 / 


1 - exp 


(12) 


A linearization  of  Eq.  (12)  and  its  substitution  into  Eq.  (7)  results 
in  an  expression  for  the  characteristic  length 

L2  = l^2  _ JS^  o . ^ (13) 

and  illustrates  how  contributions  due  to  the  set  of  consecutive  re- 
actions contract  the  penetration  depth.  This  result  is  in  agreement 
with  6rens'(14)  earlier  calculations  showing  the  effect  of  transport 
restrictions  on  the  reaction  intensity  profile.  Evidently,  the  shift 
in  the  reaction  profile  may  arise  from  any  contribution  to  the  set. 


APPLICATION  TO  Ag/AgCt  SYSTEM 

If  the  reaction  path  remains  unchanged,  the  electrode  behavior  is 
well  simulated  using  the  simple  concept  of  invariant  profile.  This  is 
always  true  for  small  time  intervals,  for  the  steady  state  or  for  the 
early  stages  of  battery  operation.  However,  in  the  course  of  battery 
operation  the  reaction  profile  is  not  invariant  either  because  of  a 
continuous  change  in  the  reaction  path  or  a change  in  the  electrode 
structure.  A description  of  this  interaction  requires  the  solution  of 
a set  of  differential  equations  whose  number  depends  upon  the  desired 
completeness  of  solution. 
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Penetration  depth-compari  son  ^^i  th_ experi ment.  In  re ce n t cotnmun i ca - 
tions(2,3),  we  examined  the  evolution  of  reaction  profile  and  found  a 
good  agreement  between  the  observed  behavior  recorded  for  the  simulated 
electrode  and  that  obtained  by  a totally  different  method  in  a structure 
resembling  actual  electrodes.  A question  thus  naturally  arises  as  to 
the  reasons  for  such  agreement.  In  an  attempt  to  answer  this  query, 
consider  the  following:  characterization  of  a porous  structure  is  a 

difficult  task  because  distribution  functions  (the  averages)  cannot  be 
determined  unambiguously.  That  is,  different  distribution  functions 
can  be  realized  for  the  same  porous  material  depending  upon  methods  and 
properties  selected  for  comparison.  On  a basis  of  rather  general  con- 
siderations, Newman (16)  showed  that  only  three  averages  appear  in  the 
treatment  of  porous  electrodes.  These  averages,  indicated  in  Figs.  2c 
and  2d,  are:  (a)  an  average  involving  the  cross  section,  i.e.,  an 

average  based  on  cutting  through  the  matrix  and  electrolyte,  <A>; 

(b)  an  average  over  the  interfacial  area,  <S>;  and  (c)  an  average  over 
the  volume,  <V>.  Thus,  similarity  in  electrode  behavior  may  be  antici- 
pated for  electrodes  having  the  same  averages. 

In  order  to  check  this  assertion  and  our  basic  premise  concerning 
the  relation  between  the  events  occurring  on  planar  electrodes  and  the 
effect  of  restraints  imposed  by  the  porous  structure,  we  have  con- 
sidered penetration  depths  experimentally  determined  for  three  vastly 
different  structures  and  by  three  different  methods  and  compared  them 
with  values  calculated  using  the  characteristic  number.  Results,  given 
in  Table  II,  show  the  anticipated  agreement  wiien  suitable  structural 
parameters  are  selected  for  systems  at  their  early  stages  of  operation. 

TABLE  II. 

Calculated  and  Observed  Penetration  Depths  for 


Various  Electrode  Structures  in  1 N KCt  at  23°C. 


structural 

penetration  depth,  cm 

structure 

parameters 

observed 

cal culated 

Segmented  electrode 

a ~ 60um 
S = 0.077  cm2 

(8  ± 0.3)10“2 

7 X 10"2 

Flotronics  elec- 
trode (membrane) 

e = 0.57 
S = 4720  cm‘^ 

(0.6  ± 0.02)10"^ 

0.55  X 10'^ 

Sphere-bed  electrode 

e = 0. 373 
S = 1010  cm~‘ 

(1.0  ± 0.2)10"^ 

0.95  X 10"^ 
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Identi fi cation  methodology.  In  this  section,  we  quote  some  of 
the  earlier  results (2,3,15)  to  demonstrate  how  the  concept  of  the  pene- 
tration depth  can  be  used  to  determine  the  probable  elementary 
processes  taking  place  at  the  onset  of  the  electrode  charging  and  how 
they  change  on  prolonged  operation. 

The  evolution  of  reaction  profile  during  the  first  minutes  of  the 
passage  of  an  anodic  current,  i(0,t)  = 50  mA  cm"^,  through  a packed 
spheres  electrode  analog,  is  shown  in  Fig.  3.  The  transfer  current 
density  distribution,  j(x,l),  is  given  in  Fig.  3a  for  one  minute  after 
the  commencement  of  electrode  operation,  while  Figs.  3b,  3c,  and  3d 
display  profiles  attained  as  indicated,  i.e.,  j(x,2),  j(x,5.5),  and 
j(x,7.5),  respectively.  In  addition,  to  accentuate  the  evolution  of 
profiles  for  diagnostic  purposes,  we  have  shown  also  two  successive 
distributions  in  Figs.  3b,  3c,  and  3d. 

In  order  to  arrive  at  the  probable  set  of  events,  consider  the 
evolution  of  the  current  density  profiles  in  terms  of  Eq.  (13).  More- 
over, assume  that  no  change  in  the  electrode  structure  has  occurred 
within  the  first  two  minutes,  a reasonable  assumption  in  view  of  Figs. 
4a  and  4b.  It  follows  therefore  that  either  the  electrolyte  conduct- 
ivity has  changed,  i.e.,  o(x,l)  > a(x,2),  or  a change  has  occurred  in 
the  dominant  step  of  the  overall  process,  thus  shifting  the  reaction 
zone  accordingly  (cf.  dashed  and  solid  lines  in  Fig.  3b).  At  a later 
time,  e.g.,  after  5 minutes  of  electrode  operation,  we  noted  the 
development  of  a "hump"  which  moves  deeper  into  the  electrode  structure 
reaching  a halfway  point  at  t = 7.5  minutes. 

The  examination  of  surface  morphology  in  the  frontal  region  within 
the  penetration  distance  of  ca  0.1  mm.  Figs.  4c  and  4d  suggests  that 
the  formation  of  a "hump"  is  due  to  blockage  of  the  second  kind(3). 

As  the  electrode  charging  continues,  the  frontal  region  becomes  heavily 
coated.  Fig.  4d,  which,  together  with  the  observed  steep  profile.  Fig. 
3d,  indicates  reaction  control  by  the  transport  of  reactant  from  the 
bulk  electrolyte. 

In  an  attempt  to  provide  additional  evidence,  we  plotted  reaction 
profile  at  comparable  times  for  a series  of  charging  currents,  namely 
for  i(0,t)  = 200,  300,  and  400pA.  The  results,  shown  in  Figs.  5a,  5b, 
and  5c,  are  in  agreement  with  the  above  analysis.  In  addition,  how- 
ever, they  point  to  the  importance  of  blocking  of  the  second  kind  for 
open  structures  and  to  the  importance  of  blocking  of  the  first  kind(3) 
for  the  electrodes  characterized  by  a large  surface  to  volume  ratio. 


CONCLUDING  REMARKS 

The  inaccessibility  to  direct  observation  usually  associated  with 
the  investigation  of  porous  electrodes  can  be  substantially  overcome 
by  the  use  of  structural  analogs.  The  segmented  electrode  arrangement 
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appears  to  be  a convenient  tool  for  following  the  evolution  of 
reaction  profiles. 

A characteristic  number,  the  penetration  depth,  can  be  formulated 
as  a product  of  kinetic  and  structural  parameters.  Such  formulation  is 
essential  for  the  interpretation  of  the  probable  operating  processes 
and  for  qualitative  assessment.  However,  additional  evidence,  e.g., 

SEM  examination  of  surface  morphology  is  desirable  to  differentiate 
between  alternatives. 

The  limitations  of  structural  analogs  as  diagnostic  tools  are 
primarily  due  to  size  of  individual  segments  and  of  spacing  containing 
the  electrolyte.  These  difficulties  appear  to  be  alleviated  by 
employment  of  methods  and  procedures  common  to  the  microelectronic 
industry. 
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Elementary  processes  for  the  case:  sparingly  soluble  re- 

actant (phase  1 1 1 )-coriducti  ve  matrix  (phase  3),  operating 
under  dissolution-diffusion  inode. 

a.  - single  pore  analog. 

b.  - set  of  consecutive  elementary  process^  . 
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Figure  2.  Porous  structure  analogs. 

a.  - real  electrode  structure  (schematic) 

b.  - electric  circuit  analog. 

' c.  - packed  sphere  analog. 

d.  - segmented  electrode  analog. 
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Figure  3.  Evolution  of  transfer  current  density  profiles.  Structure 
packed  spheres;  electrolyte  1 N KCe,  superficial  current 
density  50  mA  cm'^.  Charging  times:  a - 1 .0  minute, 

b - 2.0  minutes,  c - 5.5  minutes,  d - 7.5  minutes. 


Figure  4.  Evolution  of  surface  morphology.  Photographs  taken  at 
a - 1 .0  minute,  b - 2.0  minutes,  c - 5.5  minutes, 
d - 7.5  minutes  (cf.  Fig.  3 for  j(x)). 
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Figure  5.  Effect  of  charging  current. 

a - 200  v«A.  b - 300  pA.  c - AOO  pA. 

(Profiles  plotted  for  times  indicated.) 
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